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ABSTRACT 

We have discovered five millisecond pulsars (MSPs) in a survey of 14 unidentified FermirLAT sources 
in the southern sky using the Parkes radio telescope. PSRs J0101-6422, J1514-4946, and J1902-5105 
reside in binaries, while PSRs J1658-5324 and J1747-4036 are isolated. Using an ephemeris derived 
from timing observations of PSR J0101-6422 (P=2.57ms, DM=12pccm“ 3 ), we have detected 7-ray 
pulsations and measured its proper motion. Its 7-ray spectrum (a power law of T = 0.9 with a cutoff 
at 1.6 GeV) and efficiency are typical of other MSPs, but its radio and 7-ray light curves challenge 
simple geometric models of emission. The high success rate of this survey — enabled by selecting 7- 
ray sources based on their detailed spectral characteristics — and other similarly successful searches 
indicate that a substantial fraction of the local population of MSPs may soon be known. 

Subject headings : gamma rays: general — pulsars: individual (PSR J0101-6422) 


1. INTRODUCTION 

The Large Area Telescope (LAT, Atwood et al. 2009) 
of the Fermi Gamma-ray Space Telescope is a pulsar 
detector par excellence. Supported by a radio cam- 
paign (Smith et al. 2008) providing ephemerides for co- 
herent folding of LAT photons, the LAT has identified 
7-ray pulsations from many normal (unrecycled) pulsars 
(e.g. Abdo et al. 2010c) and millisecond pulsars (MSPs, 
Abdo et al. 2009a). Additionally, tens of new pulsars 
have been discovered in “blind” searches for periodicity 
in the LAT data (e.g. Abdo et al. 2009b; Pletsch et al. 
2012). 

A third, indirect method of detection has been ex- 
tremely successful in discovering new pulsars, espe- 
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dally MSPs 17 . Gamma-ray pulsars — both young and 
recvcled — are stable emitters and nearly all have a char- 
acteristic spectrum dN/dE oc E~ T exp(—E/E c ) with 
T < 2 and 1 < E c /GeV < 10 (Abdo et al. 2010c), 
making unidentified, nonvariable LAT sources with such 
spectra good pulsar candidates. The fine angular resolu- 
tion and large effective area of the LAT allow a typical 
source localization of ~10', about the 1 GHz beam size of 
100 m class radio telescopes. This happy coincidence en- 
ables deep pulsation searches with a single pointing. Dis- 
covery of a pulsar and a subsequent timing campaign can 
“close the loop : ’ by providing an ephemeris with which to 
fold LAT photons and resolve v ra y pulsations. Indeed, 
since most MSPs reside in binaries, initial detection and 
characterization of the orbit at longer wavelengths is the 
only way 7-ray pulsations can be detected from these 
pulsars. 

Due to sensitivity limitations, LAT-guided surveys nat- 
urally target the nearly-isotropic population of nearby 
MSPs. The deep radio exposures afforded by the effi- 
cient target selection offer the tantalizing possibility of 
completely cataloguing the local field MSPs whose radio 
and 7-ray beams cross the Earth. This relatively com- 
plete sample will be valuable in constraining the emission 
mechanisms of MSPs and in characterizing the evolution 
of their binary progenitors. 

Towards this end, we performed a small, targeted sur- 
vey at the CSIRO Parkes telescope, complementary to a 
Parkes search described in Keith et al. (2011) that had 
different source selection criteria. 

2. TARGET SELECTION 

We began with the set of southern (Decl. < —40°; 
northern sources are visible to the more sensitive Green 
Bank Telescope), nonvariable sources in a preliminary 
version of the 1FGL catalog of 7-rav sources (Abdo et al. 


17 See https:/ /confluence.slac.otanford.edu/display/GLAMCOG/Public+List- 
Detected+Gamma^Ray-t-Pulsars for an up-to-date list of LAT- 
detected pulsars 
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2010b), We further restricted consideration to sources 
whose LAT position estimates had a 95% confidence 
error radius <7', the HWHM of the 1.4 GHz beam of 
the Parkes telescope, and which had no plausible blazar 
counterpart. We classified the remaining sources by vi- 
sual inspection of the 7-ray spectrum. Sources with spec- 
tral shapes resembling known pulsars (see above) were se- 
lected for observation, while those sources best described 
by power laws or with significant spectral breaks but 
T > 2 were deemed likely to be blazars and were not 
selected. Due to the difficulty of spectral modeling and 
localization of LAT sources in the Galactic plane, we dis- 
carded sources with \b\ < 5°. Our final list comprised 14 
good candidates. The source positions (see Table 1) were 
computed from the same analysis used to estimate the 
spectral shape, but in all cases the difference in position 
with 1FGL is well within both the 1FGL 95% error con- 
tour and the Parkes beam. 

3. RADIO SEARCHES 

We observed the 14 selected sources (§2) at the Parkes 
64-m radio telescope between 2009 November 25 and 
December 8. Each source was observed for between 
1 and 2 hours (see Table 1) using the center beam 
of the multibeam receiver. We recorded 1-bit-digitized 
total-power samples every 125 /is from each of 512 fre- 
quency channels spanning 256 MHz of band centered on 
1390 MHz, writing the data to disk for off-line analysis 
(see Manchester et al. (2001) for more details on the re- 
ceiver and data-acquisition system.) 

The data were analyzed on the koala computer 
cluster at Columbia University using standard pulsar 
search techniques implemented in the PRESTO package 
(Ransom 2001). We dedispersed each data set ideally up 
to a dispersion measure of DM = 270pccm -3 , and when 
searching for periodic signals we allowed for accelerated 
signals caused by pulsar motion in a binary system. This 
was parameterized within PRESTO by the parameter 
zmax = 50, which is related to the maximum number 
of bins a signal can drift within the Fourier spectrum 
and still be properly detected by the search method. All 
of the data were analyzed in this manner within about a 
week of collection, and on 2009 December 4 we found our 
first MSP. Within a few more days we confirmed it and 
four other new MSPs. A 6th MSP was detected unbias- 
edly in our sample, but it had been previously discovered 
by Keith et al. (2011, see Table 1). 

In a few cases, the pulsars were not confirmed on the 
first attempt and required subsequent observations, ow- 
ing to their faintness and scintillation in the interstellar 
medium. We reanalyzed the data by increasing the ac- 
celeration search space to zmax = 200 and by analyzing 
only the first 30 minutes of each data file. The latter 
approach improves our sensitivity to highly accelerated 
pulsars (e.g., for an MSP in a 1-day binary, a 2-hour ob- 
servation cannot be corrected ideally under the assump- 
tion of a constant acceleration, no matter how large the 
zmax used), but no new pulsars were detected. 

The nominal sensitivity of our searches, for a 2 hr ob- 
servation and assuming a pulsar duty cycle of 20%, is 
0.11 mJy for P >2ms and DM < 40pccm“ 3 , degrading 
gradually for shorter periods and larger DMs. While in 
detail this depends on the sky background temperature, 
the contribution from the Galaxy at our observing fre- 


quency only differs by a maximum of 3 K for different 
locations, out of a total system equivalent temperature 
of nearly 30 K, so that we provide an average flux density 
limit. Two of the non-detections have limits of 0.12mJy 
and 0.14mJy owing to shorter integration times (see Ta- 
ble 1). Because of scintillation, a pulsar with a larger 
average flux than our nominal limits might not be de- 
tectable in a particular observation, and vice versa. In 
any case, the nominal flux density represents the most 
sensitive searches for MSPs done at Parkes, especially 
considering that these LAT-selected targets are largely 
expected to be relatively nearby. 

Characterization of four of the MSPs requires a longer 
radio timing campaign and/or LAT dataset than avail- 
able for this work. We defer discussion of these MSPs 
(and of scintillation effects on sensitivity, which like- 
wise requires repeated observations) to Camiio et al. (in 
preparation). The remaining pulsar presents a peculiar 
light curve and we describe it further below. 

4. PSR J0101-6422 
4.1. Timing Solution 

After the discovery observations on 2009 November 25 
and confirmation on December 8, we began regular tim- 
ing observations of PSR J01 01-6422 at Parkes. We ob- 
served the pulsar on 35 days through 2011 November 9, 
detecting it on 28 occasions; on the remaining 7 days it 
was too faint to detect, due to interstellar scintillation 
(DM=12pccm -3 .) Each observation was typically one 
hour, with the same receiver and data acquisition system 
used in the search observations. Using TEMP 02 13 with 
the 28 times-of-arrival, we obtained a phase-connected 
timing solution whose parameters are given in Table 2. 

The Shklovskii effect (Shklovskii 1970) increases the 
measured value of P = 5.2 x 10 “ 21 beyond that intrin- 
sic to the pulsar by A P = Pv\jdc , with v± the source 
velocity transverse to the line-of-sight. For nearby ob- 
jects, it may dominate the true spindown rate and lead to 
overestimates of the spindown luminosity (Camiio et al. 

1994) . The observed proper motion, at the 0.55 kpc DM 
distance (Cordes & Lazio 2002), implies a true P{ = 
(4.3 ± 0.1) x 10“ 21 . At this distance, the indicated 
v± = 41 km s _1 is typical for an MSP (see Nice & Taylor 

1995) . 

4.2. Radio Profile and Polarimetry 

We made polarimetric and flux-calibrated observations 
of PSR J0101+3422, using the PDFB3 digital filterbank. 
The data v/ere analyzed with PSRCHIVE (Hotan et al. 
2004). The best full-Stokes profile that we obtained is 
shown in Figure 1 (the flux density of this profile is not 
necessarily representative of the average pulsar intensity, 
which varies greatly owing to interstellar scintillation). 
The main pulse is linearly polarized at the ~15% level, 
but the low signal-to-noise ratio prevents a useful de- 
termination of rotation measure (the best-fit value is 
~ +10 rad m -2 , but within the uncertainties is consis- 
tent with zero). These data show that the subsidiary 
pulse component is composed of two outer peaks joined 
by a low bridge of radio emission. 

1 8 http:/ /www. atnf.csiro.au/resegtrch / pulsar / tempo2 
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Table 1 

Results of Radio Searches of 14 1FGL Sources at Parkes 


1FGL Name 

R A a 
(J2000.0) 

Decl.* 

(J2000.0) 

1 

(deg) 

b 

(deg) 

r 

E c 

(GeV) 

Obs. Time 
(hr) 

Period 

(ms) 

DM 

(pc cm" 3 ) 

Distance 11 

(kpc) 

BinaJ*y 

J0101 .0-6423 

01**00“58H 

— 64°24 , 03" 

301.2 

-52.7 

1.3 

2.3 

1.0 

2.57 

11.93 

0.55 

yes 

J0603.C-4012 

06 h 03 m 04?9 

— 40°11'02" 

246.8 

-25.9 

1.9 

6.5 

2.0 

— 

— 

— ’ 

— 

J0933.S-6228 

09 h 33 m 5S55 

-62°27'54" 

282.2 

-7.8 

0.1 

1.2 

2.0 

— 

— 

— 

■ — 

J 1036.2-6719 

10 h 36 n '16!9 

— 67°20'34" 

290.4 

-7.8 

1.5 

2.2 

2.0 

— 

— 

— 

— 

J 1227.9-4852 

I2 h 27 m 50*5 

-48° 51' 54" 

298.9 

13.8 

1.8 

1.6 

2.0 

— 


— 

— 

J 1232.2— 5118 

12 h 31 m 49?l 

-51° 18' 50" 

299.8 

11.4 

1.8 

2.3 

2.0 

— 

— 


— 

J 1514. 1-4945 

15 h 14“>05?7 

— 49°45'32" 

325.2 

6.8 

1.7 

5.9 

2.0 

3.59 

31.5 

0.9 

yes 

J1 624.0-4041 

16 h 24 m 06?2 

—40° 40' 48" 

340.6 

6.2 

2.1 

3.3 

2.0 

— 


— 

— 

J1 658.8-5317 

16 h 58 m 43“2 

— 53°17'45" 

335.0 

-6.6 

2.1 

1.8 

1.3 

2.44 

30.9 

0.9 

no 

J 1743. 8— 7620 

17 h 43 m 44“6 

-76°20'42" 

317.1 

-22.5 

1.2 

2.0 

1.5 

. — 


— 

— 

J1 747.4-4035 

17 h 47 m 29U 

-40' 36' 07" 

350.2 

-6.4 

1.5 

3.4 

1.4 

1.65 

152.9 

3.3 

no 

J1902.0— 5110 

19 h 02 m 0K5 

-51°09'43" 

345.6 

-22.4 

1.7 

4.4 

1.2 

1.74 

36.3 

1.2 

yes 

J2039.4-5621 

20 h 39 m 30?5 

—56° 20' 42" 

341.2 

-37.1 

1.6 

2.7 

1,2 


■ — 

— 

— 

J2241.9-5236 c 

22 h 41 m 52»4 

— 52°37'37" 

337.4 

-54.9 

1.6 

3.6 

2.0 

2.19 

11 

0.5 

yes 


a Parkes telescope pointing position, 

b Distances computed from measured DM using the NE2001 model (Cordes Lazio 2002). 
c See Keith et al. (2011). 



Figure 1. Polarimetric pulse profile of PSR J0101-6422 at 
1.4 GHz, displayed with 128 phase bins, based on 6 hr of Parkes 
PDFB3 data using 256 MHz of bandwidth. The black trace corre- 
sponds to total intensity, red to linear polarization, blue to circular. 
In the upper window, the position angle of linear polarization is 
plotted for bins with linear signal-to-noise ratio > 3. The pulse 
peak is displayed with arbitrary phase, and the mean flux density 
is ~ 0.2 mjy. 


4.3. Gamma-ray Profile 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Rotation Phase 


Figure 2. Phase- aligned 7- ray (blue histogram) and radio (red 
trace) pulse profiles of PSR J0101-6422. The radio profile is 
summed from 28 timing detections, while the 7- ray profile cor- 
responds to the optimal aperture described in the main text. The 
background estimation (dashed horizontal line) is derived by count- 
ing photons in the off-pulse phase window (§4.4) and is consistent 
with the off-pulse background level obtained via spectral analysis. 
The vertical dashed line indicates the phase (<j> = 0.5) from which 
the offset of “PI” is measured. 


To characterize the 7-ray profile, we selected photons 
collected between 2008 Aug 4 and 2011 Aug 1 with 
reconstructed energies 0.3-3 GeV lying within 1.1° of 
the timing position. This selection balances complete- 
ness (including as many pulsed photons as possible) and 
pulsed signal-to-noise ratio as measured by the H-test 
(de Jager et al. 1989). We apply the same data process- 
ing as described in §4.4, omitting the horizon cut to in- 
crease the livetime by about 20%, at the expense of a 
slightly increased background. 

The light curve corresponding to this extraction ap- 
pears in Figure 2. If 7 rays are produced above the null 
charge surface (the locus of points where the magnetic 
field is orthogonal to the pulsar spin axis), they appear 
to an observer in the opposite hemisphere, whereas ra- 
dio emission from low altitudes is beamed into the same 


hemisphere (see Romani & Yadigaroglu 1995). In this 
scenario, the line-of-sight to PSR J0101-6422 would in- 
tercept a cone of radio emission in one hemisphere to 
produce the weaker radio peak at <t> ~ 0.5 and a cone of 
radio and 7 emission in the other to produce the two 7 
peaks and bright radio peak at <j> ~ 1. In keeping with 
this picture and the convention established in the litera- 
ture, we identify the 7 peak at 4> ~ 0.8 as “PI” and that 
at <£~0.1 as “P2”. 

The relative phasing of the peaks is of interest (e.g. 
Watters et al. 2009). The radio-to-7 offset, 8 , is mea- 
sured from the center of the leading radio peak (taken as 

= 0.5, see Figure 2). We estimated the phase of the 
7-ray peaks using unbinned maximum likelihood with 
a two-sided Lorentzian model for the peaks. These esti- 
mates yield 8 = 0.33±0.01 8 tat±0.03 sys t and a 7-ray peak- 
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Figure 3. The best-fit model light curves for the two-pole caus- 
tic (TPC) and outer gap (OG) models (see §5.1.) The best fit 
geometries (a, C) for TPC and OG are (26°, 79°) and (90°, 36°). 

to-peak separation of A — 0.26 d= 0.03 sta t ± 0.02 syst . The 
DM-induced uncertainty in 5 is < 0.00 IF. The system- 
atic uncertainty for both S and A is estimated from the 
scatter obtained by fitting alternative functional forms 
for the peaks. 

4.4. Gamma-ray Spectrum 

To measure the spectrum of psr, we analyzed low back- 
ground DIFFUSE class “Pass 6” (Atwood et ah 2009) LAT 
data collected between 2008 Aug 4 and 2011 Apr 17. We 
filtered periods where the observatory’s rocking angle ex- 
ceeded 52° or the Earth’s limb impinged upon the field 
of view (horizon cut, requires zenith angle < 100°). We 
modelled the sensitivity to these events with the flight- 
corrected P6_V 11 .DIFFUSE instrument response function. 
To compute background contributions, we used a prelim- 
inary version of the 2FGL catalog and its accompany- 
ing diffuse emission models 19 . We selected events within 
10° and re-fit the spectra of point sources within 8° and 
the normalizations of the diffuse models. Fits employed 
pointlike (Kerr 2011), a binned maximum likelihood al- 
gorithm. 

The phase-averaged spectrum of PSR J01 01-6422 is 
well described by an exponentially cutoff power law, 
dN/dE = No (E/Eo)~ r exp(— E/E c )^ whose parame- 
ters appear in Table 2. We also considered separately 
the spectra of the three major components, viz. P2 
(0 < (j> < 0.33), the off-pulse (0.33 < (j) < 0.63), and 
PI (0.63 < 0 < 1). The parameter uncertainties are 
relatively large, and accordingly we found no significant 
difference in the spectral shapes of the two peaks. Fur- 
ther, we detected no significant emission in the off-pulse 
phase v/indow and derived a 95% confidence upper limit 
on the energy flux of a point source with a power law 
spectrum with T = 2.2, obtaining 6 x 10 _12 erg cm -2 
s _1 (>100 MeV; scaled up by 3.3 to the full phase win- 
dow.) 

The observed flux is related to the 7-ray luminosity 
by Fj — L 7 /47r/nd 2 , where fa is an unknown beaming 
factor ' typically of order unity, e.g. Watters et al. 2009) 
implying a phase-averaged, isotropic (fa = 1) 7-ray lu- 
minosity fa = 4.0 x 10 32 ergs _1 = 0.04 E at d = 0.55 

19 http : / /fermi .gsfc.nasa.gov/ssc / data/access/lat /Backgrounds 


kpc. Best-fit values of fa from the TPC and OG models 
(§5.1) are 0.7 and 0.9. This 7-ray efficiency is in accord 
with other LAT-detected MSPs (Abdo et al. 2009a). 

4.5. X-ray and Optical Observations 

A 3.3 ks Swift observation with the X-ray Telescope in 
PC mode was obtained on 2009 Nov 21 to search for an 
X-ray counterpart to the then-unidentified 1FGL source. 
No significant counterpart is detected, with a 3 a upper 
limit at the pulsar position of 2.0 x 10 -3 ctss“ 1 (0.5-8 
keV). Assuming a power-law spectrum with photon index 
T = 1.5 for a column with iVn = 3.7 x 10 2O cm -2 (10 
HI atoms per free electron along the line-of-sight) yields 
an unabsorbed flux limit of 1.1 x 10~ 13 erg cm^ s' 1 , or 
efficiency Lx/E < 5.0 x 10 -4 at 0.55 kpc. This limit is 
comparable to the detected X-ray flux from other MSPs 
(e.g. Marelli et al. 2011). 

We searched four DSS plates and found no optical 
counterpart consistent with the pulsar’s position, sug- 
gesting m > 21. 

5. DISCUSSION 
5.1. PSR J0101-6422 

At first blush, PSR J0101-6422 is an unremarkable 
member of the growing population of 7-ray MSPs. The 
orbital characteristics suggest an epoch of Roche lobe 
overflov/ from an evolved 1-2 M© companion which later 
became the M).2M© He white dwarf currently inferred 
in the system (e.g. Tauris & Savonije 1999). However, as 
we discuss below, PSR J0101-6422’s light curve is chal- 
lenging to explain with simple geometric models that tie 
the radio and 7 emission to particular regions of the mag- 
netosphere. 

Remarkably, a simple picture of 7-ray emission aris- 
ing in the outer magnetosphere appears to describe the 
majority of both young (though see Romani et al. 2011) 
and recycled pulsars (e.g. Venter et al. 2009). On the 
other hand, while radio emission in young pulsars and 
some MSPs appears to come from lower altitudes, sev- 
eral MSPs have phase-aligned 7 and radio peaks (e.g. 
Abdo et al. 2010a), implying a joint emission site at high 
altitude. These observations raise the prospect of MSPs 
whose radio emission combines “traditional” polar cap 
emission with a higher-altitude component. 

To determine which scenario best describes the emis- 
sion pattern of PSR JO 101-6422, we performed joint 
fits to the radio and 7-ray data following the method 
of Johnson (2011). We modelled the radio as a sim- 
ple ‘ cone” component (Story et al. 2007) and the 7 
rays using both the “two-pole caustic” (TPC) model 
of Dvks & Rudak (2003, N.B. we adopt a larger max- 
imum cylindrical radius of 0.95 of the light cylinder) 
and the “outer gap” (OG) model (Cheng et al. 1986; 
Romani & Yadigaroglu 1995). The angle between the 
pulsar spin axis and the line-of-sight, £, is unknown a 
priori , but if the radio emission arises from low altitudes, 
then the presence of two peaks separated by roughly half 
of a rotation indicates the inclination of the magnetic 
axis from the spin axis (a) cannot be small. The best-fit 
light curves for the two 7-ray models appear in Figure 3. 
The low a = 26° preferred for the TPC model cannot 
produce the radio interpulse and produces too much off- 
j^Sl&e 7-ray emission. The OG model, with a best fit at 
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Table 2 

Measured and Derived Parameters for PSR J0101-6422 


Parameter Value 

Right ascension, R.A. (J2000.0) 01 h 01 m 11® 1163(2) 

Declination, decl. (J2000.0) -64°22'307 171(2) 

Proper motion in R.A. X cos decl. (mas yr -1 ) .... 10(1) 

Proper motion in decl. (mas yr -1 ) —12(2) 

Galactic longitude, l (deg.) 301.19 

Galactic latitude, b (deg.) -52.72 

Spin period, P (ms) 2.5731519721683(2) 

Period derivative (apparent), P 5.16(3) X 10 -21 

Epoch of period (MJD) 55520 

Orbital period, P b (days) 1.787596706(2) 

Epoch of periastron (MJD) 55162.4011764(3) 

Semi-major axis (It s) 1.701046(2) 

Eccentricity < 1.5 X 10“ 5 

Dispersion measure, DM (pc cm -3 ) 11.926(1) 

Spin-down luminosity 0 ”, E (ergs -1 ) 1.0 X 10 34 -16% 

Characteristic age a , r c (yr) 9 X 10 9 +20% 

Surface dipole magnetic field strength (Gauss) a ... 1.1 x 10 7 -9% 

Mean flux density b at 1.4 GHz, Si,4 (mJy) 0.28 + 0.06 

Radio-^y-ray profile offset, 5 (P) 0.35 ± 0.01 ± 0.05 

7-ray profile peak-to-peak separation, A (P) 0.26 ± 0.03 ± 0.02 

7-ray (> 0.2 GeV) photon index, T 0.9 ± 0.3 + 0.2 

7-ray cut-off energy, E c (GeV) 1.6 + 0.4 + 0.3 

Photon flux (> 0.1 GeV) (10 -9 cm -2 s -1 ) 9.5 ± l-fij^o 

Energy flux (> 0.1 GeV), F 7 (lO"” 11 erg cm -2 s _1 ) 1.1 ±0.1 ±0.1 


Note. — Timing solution parameters given in TDB relative to the DE405 plane- 
tary ephemeris. Numbers in parentheses represent the measured 1 a TEMP 02 timing 
uncertainties on the last digits quoted. For 7-ra” parameters, the first uncertainty is 
statistical and the second systematic. 

a Corrected for Shklovskii effect at 0.55 kpc; the percentage change from the nominal 
value is given. 

b Average flux for 28 timing detections. 


a = 90° (orthogonal rotator), produces two radio peaks 
but fails to produce two clear 7 peaks and the radio in- 
terpulse morphology. 

Since neither model faithfully reproduces the observed 
light carves, we propose that one or both of the radio 
peaks may originate at higher altitude. Emission over a 
range of altitudes would also explain the low observed 
linear polarization of PSR J0101-6422 (e.g. Dyks et al. 
2004) and admit smaller a values. Alternatively, the fail- 
ure of geometric models may indicate that the physical 
details of the MSP magnetosphere—e.g., currents, mul- 
tipolar fields, and plasma loading of the field lines — and 
the radio and 7 emission mechanisms play an important 
role. 

5.2. Survey Implications 

With 6/14 pointings resulting in unbiased MSP detec- 
tions, this survey was fruitful, though this high efficiency 
may be luck. Indeed, in our timing campaign, we only de- 
tected PSR J0101-6422 on 28 of 35 attempts, and several 
observations required prior knowledge of the ephemeris 
for detection. Thus, because of scintillation, we estimate 
that PSR J0101-6422 is only detectable in 2/3 observa- 
tions, depending somewhat on integration time. Many 
field MSPs are in eclipsing binaries, further diminishing 
the probability of detection in a single observation. 

We believe the strongest factor in our search’s success 
was the source selection criteria. E.g., Ransom et al 
(2011) selected nonvariable LAT sources (but includ- 
ing no spectral shape information) and found MSPs in 
3/25 pointings. But when considering only targets with 


\b\ > 5°, their efficiency jumps to 3/8. Keith et al. 
(2011) employed some spectral information and detected 
MSPs unbiasedly in 1/11 pointings (1/4 for |6| > 5°). 
A short list of 7 spectrally-exceptional candidates drawn 
up for a Nan^ay radio telescope survey included 5 MSPs, 
three of which were discovered from Nanqay observa- 
tions (Cognard et al. 2011; Guillemot et al. 2011). And, 
Hessels et al. (2011) used nearly identical selection crite- 
ria as this work in a Green Bank Telescope survey at 350 
MHz and found MSPs in 13/49 pointings, all but one at 
|6| > 5°. the 

Targeting pulsar- like sources is clearly efficient. How- 
ever, new LAT sources will lie at the sensitivity threshold 
and the limited statistics will not admit classification of 
spectral shape and variability. Nonetheless, we argue 
that future radio searches should omit detailed ranking 
schemes and target any LAT source off the Galactic plane 
that has not been associated with a known blazar. 

We believe this substantially increased target list (and 
telescope time) is justified. In the 2FGL source catalog 
(Abdo et al. 2011), there are ~350 unassociated sources 
more than 5° off the Galactic plane. MSPs account for 
roughly 5% of associated LAT sources, so 15-20 MSPs 
may remain to be found in these unassociated 2FGL 
sources. With only 350 positions to monitor, multiple 
deep observations can ameliorate the confounding effects 
of scintillation and eclipses while requiring a fraction of 
the time of a comparably complete all-sky survey. By 
comparison, the otherwise prodigious Parkes multibeam 
survey found “only” 17 MSPs in over 40,000 pointings 
(e.g. Faulkner et al. 2004). In addition to the intrinsic 
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interest in discovering new MSPs, by thoroughly search- 
ing all such LAT sources, we will gain confidence that the 
MSPs we have found can truly provide a volume-limited 
sample of the energetic, 7-loud MSPs, an invaluable step 
in the study of the MSP population of the Milky Way. 
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